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Brezin-Gross-Witten model as "pure gauge" limit of 

Selberg integrals 

A.Mironov* A.Morozov^ and Sh.Shakirov* 
Abstract 

The AGT relation identifies the Nekrasov functions for various AT = 2 SUSY gauge theories 
with the 2d conformal blocks, which possess explicit Dotsenko-Fateev matrix model (/3-ensemble) 
representations the latter being polylinear combinations of Selberg integrals. The "pure gauge" 
limit of these matrix models is, however, a non-trivial multiscaling large- N limit, which requires 
a separate investigation. We show that in this pure gauge limit the Selberg integrals turn into 
averages in a Brezin-Gross-Witten (BGW) model. Thus, the Nekrasov function for pure SU(2) 
theory acquires a form very much reminiscent of the AMM decomposition formula for some model 
X into a pair of the BGW models. At the same time, X, which still has to be found, is the pure 
gauge limit of the elliptic Selberg integral. Presumably, it is again a BGW model, only in the 
Dijkgraaf-Vafa double cut phase. 

1 Introduction 

CD ' 

r"" 1 1 ' The pure gauge limit. The AGT relation [l]-[50] is an explicit formulation of duality between 
the 2d and 4d descriptions of conformal 6d theory of self-dual 2-forms, compactified on a Riemann 
surface [52]. The theory of the corresponding M5-brane is long known to be related to integrability 
theory [53], but explicit route from integrability to the AGT relation still remains a mystery. 

A promising approach to origins of the AGT relations is through their reformulation as relations 
■ between matrix models [54] and Seiberg-Witten theory [55, 56], see [57] for a concise review of this 
idea (which is a new application of the topological recursion [58]-[61]). 
t-h ' Surprisingly or not, despite obvious conceptual advantages of such an approach, some simple 

properties of original AGT relations are not so easy to describe in the matrix-model reformulation. 
A typical example is the "pure gauge" limit (PGL), where the dimensional transmutation takes place 
and the conformal invariance gets broken. In the matrix model formulation, this corresponds to a 
^ non-trivial double-scaling large N limit of the relevant matrix models, which will be the subject of 
the present paper. 

In this paper, we concentrate on the simple case of pure SU(2) Nekrasov function, which 
(in the context of the AGT relations) arises as the PGL of either the 4-point conformal block 
_B(°)(Ai, A2, A3, A4; A, c\x) on sphere [7] or the 1-point conformal block (A ext ; A, c\q) on torus 
(where q = e 2mr and r is the modulus) [21, 51]: 



(N 
> 



£*(A|A) = lim r J B(°)(A 1 ,A 2 ,A3,A 4 ;A,c|x)l = lim \ B& (A ext ; A, c\q) 



ID 



a;(A2-Ai)(A3— A 4 )=A 4 ?A^ t =A4 
Our aim is to study this pure gauge limit at the level of matrix models. 
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PGL in the matrix model formulation. Fortunately, matrix model (i.e. the Dotsenko-Fateev- 
like /3-ensemble) representations are already known both for B^ [12, 51] and for [46]. The first 
one is represented [36] as an AMM decomposition [59, 61] into two spherical Selberg integrals 



*«»(*) =exp [-- p j: t 



2 ^ x k \v + d 



2 + dti 



v- d 



k J 



zP(u + ,v + ,N + \t + )z£>(u_,v_,N^\t-) (2) 



'(0), 



while the second one is the single, but elliptic (toric) Selberg integral 



BW (q) = f[(l - qT (A, N\q) 
p=i 

Then, eqs.(l), (2) and (3) imply that 

5*(A|A) = exp ("Ig^^y ^(* + )^(*")L = Z ^ 



(3) 



(4) 



and this is the straightforward matrix model formulation of eq.(l). The partition functions Z^_ 

and denote here what we call the PGL of the corresponding Selberg integrals. To describe this 
PGL, explicit expressions for the Selberg integrals are needed. Note that the second equality in (4) 
automatically provides an AMM decomposition formula for Z^ into a pair of the zf^ models, even 
before explicit expressions are discussed. 



The Selberg integrals. Explicitly, the Selberg integrals have a form of the eigenvalue /^-ensembles 



N N 

Zf (u, v,N\t) = ± fd Zl ... dz N H( Zl - z 3 r ]I " 1) *«f 



(5) 



2?r N N 

zg ) (A,N\q)=j-[ dz 1 ...dz N f[6*(z i -z j \q) 2f3 f[e*(z i \q) 



(6) 



where / = \f—l and 9*(z\q) = sin(z/2) — qsm(3z/2) + . . . is the normalized odd theta-function on 
torus [46] . The normalization constants Sq , Si are needed to satisfy the requirements zf^ (t = 0) = 
Z < g\q = 0) = 1 implied by the conditions B^°\x = 0) = B^\q = 0) = 1 for properly normalized 
conformal blocks. It is these /3-ensembles that we will use to study the pure gauge limit. 

Parameters of Selberg integrals. To study the PGL in terms of the Selberg integrals, one needs 
to describe clearly the values of their parameters. In the spherical case [12, 51], they are given by 



a 1 -a 2 + a-e 1 -e 2 ai - ei - e 2 n a 2 
N + = , u + = 2 , v + = -2— 



(7) 



«4 — «3 — a 



u_ = 2 



«4 — ei — £2 

£2 



w_ = -2 



«3 
«2 



(8) 
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and, in the toric case [46], they are given by 



at aext . 2a + ei + e 2 _ 2a 

7V = , A= = — , u = 3A ext + 3N-l (9) 

where the a-parameters are related to the initial A-parameters (conformal dimensions) via 

A(a) = " (6l + 62 " a) (10) 
ei£2 

and 

c= i_ 6 ffi±^y ( n) 

is the central charge. Here it is convenient to write all these formulas in terms of the Nekrasov 
parameters ei, e 2 , which are in one-to-one correspondence with the matrix model parameters j3 (the 
power of the Vandermonde determinants) and g s = g (the "string" coupling constant, aka the genus 
expansion parameter): 

ei = -gy/p, e 2 =g/y^ (12) 

or vice versa 

g 2 = -eie 2 , (3 = -^ (13) 

62 

This completes the list of relations between the parameters, and allows one to look at the PGL of 
the Selberg integrals. This limit is simple in terms of the external dimensions: 

ai — > 00, afx = fixed = A 4 (14) 

a e xt — > 00, a^ xt q = fixed = A 4 

while in terms of the matrix model parameters it gets more sophisticated. We will now describe the 
limit in terms of N, u, v (for the spherical Selberg models) and of N, A (for the elliptic model). 

PGL of Selberg integrals. Relations (7)-(8) imply that, in the PGL, the parameters u, v, N 
of the spherical Selberg integrals all tend to infinity. However, the same relations indicate that a 
particular combination of parameters, that is, u + v + 2/3N remains finite in the PGL, since it does 
not depend on the external dimensions. We find it most convenient to parametrize this combination 
by a single variable n 



which is equal to 



PGL(w + v + 2/3 AO = /3n + /3 - 1 (15) 



_L 2a ei + e 2 

n± = ± — , a = a — (lb) 

ei 2 



for our "+" and "-" Selberg models. Consequently, the PGL for the spherical Selberg model looks 
like a non-trivial double scaling limit, where the parameters u,v,N tend to infinity as 



u+v+20N=/3n+/3-l 



Z?'(n\t)= lim Z$>(u,v,N ( ^ + ^ 2)fc ) (17) 
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where distinguished combinations (uN + /3iV 2 ) + and (uN + /3/V 2 )_ in the PGL are proportional to 
[i\[i2 and /U3/U4, respectively (/x's are mass parameters of the Nekrasov function, linearly related to 
a's). This particular rescaling of the time-variables 1— > tk(uN + f3N 2 )~ k is necessary to suppress 
a growth of correlators in the model: only with variables defined in this way, the partition function 
has a finite PGL. In particular, only with such a rescaling of variables the decomposition formula 
(2) remains non-trivial in the PGL (14) and, moreover, turns into formula (4), i.e. 

Zi°2(t±) = zi°\n ± \t ± ) (18) 

Similarly, in the toric case, relations (9) imply that the PGL for the toric Selberg integral is 

zi 1) (A\A)= lim ZP(A,N\q) (19) 

N— too, q— >0 

q/3 2 iV 4 =A 4 

where, since no time- variables are introduced, no additional rescalings are required. As one can see, 
the PGL's of the Selberg models are quite sophisticated: it is by no means transparent that eqs.(17) 
or (19) do at all have a finite limit. However, as we shall see below, they do, and the main problem 
is to give some constructive description of this limit. This paper is devoted to finding a (at least, 
partial) solution to this problem. 




PGL of spherical Selberg: BGW model. As the first (simplest) part of solution to this 
problem, in this paper we demonstrate that zi°\ the PGL of the Selberg partition function Z^ is 
actually the partition function of the (/3-deformed) celebrated BGW model [62, 63, 61] of size n and 
in the character phase [63] : 

/3 2 (/3n + 2-2/3) 2 ^ 

+ (/3n + 1 - /3)(/3n + 2 - /3)(/3n + 1 - 2/3) ~ (/3n + 1 - /3)(/3n + 2 - /3)(/3n + 1 - 2/3) * 2 + " ' 

where the integral over U is the /3-deformed unitary integral, and Vol^n) is the /3-deformed volume 
of the unitary group. As usual for the BGW model, the time-variables are identified with traces 
of the external field powers tj~ = tr^ k /k, and this brings us directly to the topic of /3-ensembles 
with external fields, which is somewhat underinvestigated and not exhaustively covered in the ex- 
isting (physical) literature. The point is that the /3-deformations are usually defined for integrals 
of eigenvalues only. Whereas the notion of trace (of determinant, etc.) in (20) remains well-defined 
as combinations of the eigenvalues, the treatment of the external field term trU^ in (20) deserves 
some comments. Actually one needs only the [/-integrals (averages) of such quantities, they will be 
defined in s.2 with the help of a /3-ensemble version of the Harish-Chandra-Itzykson-Zuber integral. 

In this paper we prove eq. (20) in two independent, but complementary ways. After the /3- 
unitary BGW model is defined in s.2, in s.3 we demonstrate that its Jack expansion (the /3-ensemble 
counterpart of the character expansion) coincides with the PGL of the Jack expansion for the spher- 
ical Selberg model. This is just an algebraic exercise, which still may be not too much transparent. 
A more conceptual way may be the method of s.4, where one instead takes the PGL of the Virasoro 
constraints (Ward identities) for the spherical Selberg model, and then shows that they coincide 
with the known Virasoro constraints [63, 61] for the character phase of the BGW /3-ensemble. 

Formula (20) fully defines the middle part of formula (4). It still remains to explain how the 
time derivatives can be taken in the external field BGW model: the simplest possibility is provided 
by a Fourier-like transform in the /3-character calculus, which is reminded in s.2. 
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PGL of elliptic Selberg: double BGW model. The second, harder part of the solution, 
taking the PGL of the elliptic Selberg model is not completely finalized in the present paper. The 
corresponding partition function Z*^(A|A) should be given by some /3-ensemble with the partition 
function 

2/3A 4 /3 2 W 2 .4 2 + /3-8(/3-l) 2 ) 

Z£\A\A) = 1+ P +- ^ > v + (21) 

A ~ ~ !) \^A 2 - (p - 1) 2 J [A 2 - (2/5 - 1) 2 J (A? -(13- 2) 2 J 

Note that the problem of finding the PGL of the elliptic Selberg model is just the same as finding 
the PGL B*(A|A) of the conformal block/Nekrasov function. Having in mind the general context 
of the problem and lessons from the (successful) solution in the spherical case, one can go further in 
several directions. 

One way to solve this problem is to attack it directly from the elliptic side, by trying to take the 
Inozemtsev limit [64] of the Ward identities for the elliptic Selberg model This remains to be 

done. 

Another way is to make a direct educated guess for what the /3-ensemble in question should be. 
Such an attempt has actually been made long ago in [65]. The conjecture was that is again a 
BGW model, but this time in another phase: the double-cut DV phase, i.e. 

(1) presumably D y 2 
Zj* — A BGW 

This suggestion can be checked, for example, by a derivation of AMM decomposition of Z^? w and 
comparison with the second equality in (4). This also remains to be done. 

The third possibility is to apply the character calculus of [66], either to the elliptic Selberg 
integral itself or to the decomposition formula (4). This is what we do in the present paper, in 
s.5 below. Applying the character calculus to the decomposition formula (4), one obtains a double 
/3-unitary-ensemble 



(23) 



where ksk = tr U k ,ksk = trU k and m± = n±+(fi— Thus, Zbgw's m the integrand are actually 
functions of U and U, while their conjugates and W enter through the mixing (intertwining) 
determinant. Note that the /3-unitary integrals in (23) have sizes n± = ±2a/e\, while the BGW 
models in the integrand have sizes m± = ±2a/e\. As explained in s.5, this shift of sizes can be seen 
as a natural property of the Fourier transformation for the /3-ensembles. 

Eq.(23) or, at least, the first terms of its A-expansion (21) can be checked in practice by substi- 
tuting (20) into (23) and taking the remaining averages 

n n n 

j f(U)[dU]f, = j #i...y # n /(e^,...,e^)n|e^-e^| 2/3 (24) 

nxn —it —it a< b 

which is valid when f(U) is an invariant function (depends only on eigenvalues or, what is the same, 
on traces of powers of U). The averages, necessary to reproduce the first two orders in A, are 

f ^ -trtftr(^) = - " , / -Mk.tr ^tr(^)- = to ~ ^ +/ ' , (25) 



Vo\p{n) y ' Pn + P-V J Volp(n) v ' [fin + P)(Pn — l)(/3n + P — 1) 
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l^trC/trC/tr(C/+) 2 = / -MkUr Uhr (t/+)tr (£/+) = '"^/j „ (26) 



Vo\ p (n) v ' J Vol^(n) v y v y (/3n + 0)(0n - l){ftn + ft - 1) 



/ 



i^(trt/) 2 (tr^) 2 = - 5 - (27) 



Vo\ fj (n) y Jy ' {ftn + ft)(ftn-l){ftn + ft-l) 



Using them, one easily verifies that (23) reproduces (21), so that operationally the model is well- 
defined and gives correct results. However, there are many conceptual questions left, in particular, 
the relation of the model (23) to the Inozemtsev limit of the elliptic Selberg model and, most 
importantly, to the double-cut BGW integrals. If the DV conjecture (22) is true, the model (23) 
should probably be interpreted as an integral representation of Z^q W . We do not consider this topic 
in the present paper, this will be done elsewhere. 

2 Comments on the definition of /^-deformed BGW model 

Generalization of unitary integrals to ft / 1 deserves comments. The problem is to give some 
concrete definition for the /3-deformed unitary integral of the form 



/ 



f(U)[dU] p (28) 



which for ft = 1 is the well-defined integral over the compact Lie group U(n) with the Haar invariant 
measure [dU\p = \. We are unaware of any similar group theory definition for generic ft (occasionally, 
such a definition exists for ft = 1/2 and ft = 2 in terms of the groups 0(n) and Sp(n), respectively). 
Instead, various other, more or less natural definitions can be suggested. 

The simplest case is when f(U) is an invariant function, i.e. it depends only on the traces of 
powers of U or, equivalently, only on the eigenvalues of U. In this case, the most natural definition 
of /^-generalization, motivated by consideration of the above-mentioned three group integrals (U(n) 
for ft = 1, 0{n) for ft = 2 and Sp(n) for ft = 1/2) is the following eigenvalue integral: 

n n n 

j f(U)[dU] p = j dfa.-.J d4> n /(diag[ e ^ 1 ,...,e 7 ^])n| e / ^-e^| 2/3 (29) 

nxn —n —n a< ^ 

i.e. the role of the /3-deformation is just to have the power 2ft of the Van-der-monde determinant. 
It is this definition which is most commonly recalled when the words " /3-ensemble" are mentioned. 

However, for the purpose of present paper this definition is not enough. The case when the 
integrand f(U) depends only on the eigenvalues, does not cover all the eigenvalue models [69, 54]. 
In particular, the main object of the present paper, which we use to describe the PGL of the Selberg 
integrals, has an integrand which is not a function of the eigenvalues of U only, it involves an 
"external field" matrix ^ in the following way 

ww= / Mk-A**«* w ) (30) 

J V0L3 (n) 

nxn 

so that definition (29) is not applicable. Note that for ft = 1 the integral Zbgw{^) obviously 
depends only on the eigenvalues of the matrix ^. We can quite naturally assume the same property 
to hold for all ft. After that, the two options still remain: to consider the integral as a function of 
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traces of positive powers ktk = tr^ k or of negative powers krt = tr^~ k . As is well-known in the 
theory of the ordinary (/3 = 1) BGW model, these two choices lead to different results, commonly 
known as the character phase ZrgWc^) and the Kontsevich phase ZsGWkiT), respectively [63]. In 
this paper we are interested in the definition of the first one, ZsGWc(t)- 

Our definition refers to the Harish-Chandra-Itzykson-Zuber (HCIZ) integral over the unitary 

matrix V, Z IZ (t,s) = f - . e ptT * vt/yt which is a function of kt k = tr^ k and ks k = tiU k . 

nxn Vol^ (n) 

To see this at /3 = 1, it is enough to diagonalize the matrices, \l/ = V^^dV-q, an d U = VijUdYy and 
use invariance of the Haar measure [dV] to change V — > V^VVjj- Then 

ziz(t,s) = f W e **vuvt = f JM ^wm* (31) 

p=i J Vol(n) J Volp(n) 

nxn nxn 

is indeed a function of the eigenvalues {^i} and {Uj} only, and thus of t and s (the symmetry under 
permutation of the eigenvalues is obvious). Once Zjz is defined, one can write 

Z B GWc(t)= [ ^Y-,Z IZ (t,s)e^ ut (32) 
J Vol/3 (n) 

nxn 

which is now a function of the positive time-variables t. Moreover, the integrand also depends only 
on the eigenvalues of the integration matrix-variable U, so that definition (29) is applicable. 

Thus, to define Zbgwc{P), it suffices to define the Itzykson-Zuber integral (31) in some indepen- 
dent way. Such a possibility is provided by the character calculus [66]: for j3 = 1, the IZ integral 
can be defined as an expansion 

Ziz(t,s) = J2j^f-,XR(t)XR(s), /3 = 1 (33) 

where \R are the characters, and d,R,Dn are their values at particular points oIr = XR^-k = $k i)j 
Dr{ti) = XRitk = n/k). The quantities oIr and Dr have an important representation theory 
meaning: they are dimensions of representations of the symmetric- and GL{n)- groups, respectively 
(labeled by the Young diagram R). Thus, for arbitrary j3, one can naturally define 

Ziz(t, s) = J2 j^3R(t)j R (s), V/3 (34) 

where jr are the well-known /3-characters (actually, the properly normalized Jack polynomials [67], 
reviewed in detail in Appendix 1 of the present paper) and dR, Dr are their values at particular 
points dR = jRifu = <5fc,i)> DR{n) = jR{tk = n/k). For arbitrary f3, these quantities do not have 
any straightforward representation theory meaning, they can be only thought of as /3-deformations 
of the dimensions of symmetric and general linear groups. Note that the n-dependence in formula 
(34) emerges only due to the n-dependent quantity Dr{ji). 

Eqs. (32) and (34) provide our constructive definition of the BGW partition function in the 
character phase. It should be emphasized, that in this way the /3-deformation of any unitary 
eigenvalue model with no more than one external field in character phase can be defined. 
This is enough for our purposes here, but not enough in principle: to study the /3-deformations of 
the Kontsevich phases or the /3-deformations of the models with multiple external fields (known 
as non-eigenvalue models), some other definitions have to be invented. Hopefully, there exists a 
unifying framework for the /3-deformations of any matrix model, not necessarily of the eigenvalue 
type. This framework remains to be discovered. 
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Note that definition (34) already appeared in mathematical literature, see, for example, [68]. It 
goes without saying that it reproduces the HCIZ integrals not only for unitary ((3 = 1), but also for 
orthogonal (/? = 1/2) and symplectic (/3 = 2) matrices. 

Given such a definition of ZbgWc /3-ensemble, it is possible to study if it indeed coincides with 
the PGL zi°^ of the spherical Selberg /3-ensemble. The next two sections are devoted to two ways 
of proving this. The first method makes a direct use of properties of the /3-characters jn, while the 
second method relies upon Ward identities for the /3-ensembles. 

3 Eq. (20) via Jack expansion 

To check the equality zi°^ = ZbgWc we expand the both quantities in the basis of Jack polynomials. 
The right hand side. For ZbgWc eqs.(32) and (34) imply that 

ZsGWcit) = J Z IZ (t,s)e^ W [dU]p = Y^jRit)^- 1 3R{s)e^ u \dU}p (35) 

To deal with the r.h.s. integral, one needs just two properties of the /3-characters: the completeness 



ex P ktk * k = JZini^jRit) 

V k J R 



(36) 



and the Haar orthogonality 

/ jR{s)j k {s-')[dU], = ^R^^y ^ = ^=± (37) 

where n is the size of U. The first identity (with ktk = tr (U^) k , ktk = Sk,i) implies that 

=^d R j R (s-') (38) 
R 

and the second property allows one to perform the integration over U : 

Z BG Wc(t) = J Z IZ (t,s)e^ u \dU]p = EEJ«WS / m^JR^ldU]? (39) 

R R 



?5^fl* w (40) 



The expansion obtained 



di 



Ft 



(41) 



is of course a /3-deformation of the well-known character expansion [66] 

Z BG Wc{t) = Y, -^n) XR{t) ' P = 1 (42) 

but not quite a naive one (because of the non-trivial 5-shift in the denominator). Note that the 
matrix sizes enter the character expansions only through the explicit factors Dr. Let us now derive 
a similar expansion for the Selberg partition function. 
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The left hand side. The Jack expansion for the Selberg /3-ensemble 

N N 



zP(u,v,N\t)=^Jdz 1 ...dz N f[{z i -z j ) 2 ^z^z i -l) v e^^ t ^ (43) 



i<3 i=l 

using the same completeness condition (36) can be written as 

,(o) 



Z ( °\u,v,N\t)=Y J {jR)3R(t) (44) 



R 

where the coefficients (if?) are averages of the Jack polynomials of the z- variables, 



JR 



1 / \ N N 

j-fd Zl ... d ZN j R u = e z */ k - ^ n - i y ( 4s ) 

n V i / i<j i=l 



Fortunately, the averages of Jack polynomials in the Selberg model are well-known to be simple 
quantities [70] . They factorize nicely, and can be generally expressed by the Kadell formula 

/ . \ = <h_ [N(3] R [u + (N -I)f3 + 1] R 

W p\ R \ [u + v + 2N(3 + 2-2(3] R 1 ' 

where notation \...\r (see Appendix 1) stands for 

= ■ ^ = II (x-/3(<-l) + (j-l)) (47) 
R (ij)efl 



Using this Kadell formula, it is immediate to take the PGL (15): 



JR 



d 



R 



* D R (n - 5) 



(uN + (3N 



2\\R\ 



(48) 



As one can see, the correlators grow like (uN + (3N 2 )\ R \ in the PGL. This is actually the reason to 
introduce the corresponding rescaling of time-variables in the PGL of the partition function: 



Zi 0) (n\t)= lim Z£'[u,v,N 

u ,v ,N — >oo 
u+v+2l3N=f3n+j3-l 

Substituting here eq.(48), one finds 



,(o) 



tk 



(uN + @N 2 ) k 



E 

R 



JR 



3R{t) 



*(uN + fiN 2 )\ R \ 



(49) 



z| 0) (n|*)=E 



4 



R 



D R (n - 5) 



3R(t) 



(50) 



what precisely coincides with (41). In this way, the PGL limit of the Kadell formula reproduces the 
BGW model: we conclude that Z^ = ZbgWc, since their Jack expansions are just the same. Let 
us now pass to the second method of proving this statement. 
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4 Eq. (20) via Virasoro constraints 



Another way to see that = Zbgw is to study the PGL of the Virasoro constraints for the Selberg 
model. Just as any matrix model (or the /3-ensemble) , the Selberg model can be characterized by 
certain linear differential equations, which arise as a consequence of the reprarametrization invariance 
of the multiple integral (i.e. as the Ward identities). In the case of the Selberg model, these linear 
differential equations have the form 



K m K+m a+b=k a h=l 



Zs(h,t2, ■ ■ ■) 



= P\N+J2ki + uN + N(N-l)pjZs(t 1 ,t 2 ,...), k>0 (51) 

Their derivation is given in Appendix 2. These equations completely determine the partition function 
of the Selberg model, therefore, one can take the PGL in these equations, not in the integral. In the 
PGL, all the three parameters u, v, N of the Selberg models tend to infinity in such a way that their 
combination u + v + 2(3N is held finite. In this limit, the Virasoro constraints get simplified (many 
terms can be thrown out) and turn into 



(„ + , + m _ {k + m _ „) A. + ,, yj mtm JL_ + g + ph 

for any k > 0, and these are precisely the Ward identities for integral (32). They can be found in 
[63, 61], of course, only for the most popular case of (3 = 1. Despite there are no doubts that eqs. 
(52) hold for arbitrary f3, in principle, it would be nice to derive them directly from the definition of 
the BGW /3-ensemble (in the present paper, the role of such a definition is played by the character 
calculus). This remains to be done. 

5 Eq. (23) via Fourier transform 

Definition (41) can be directly used to convert the action of the intertwining operator in (4) into 
a kind of a 2-matrix BGW model. This can be done via the Fourier transform, widely used in the 
character calculus [66]. The basic and most important Fourier relation has the form 

'«M-^/^*M (53) 

nxn \ k / 

where ksk = tr U k . This relation is a direct consequence of the completeness condition (36) and the 
Haar orthogonality (37): one substitutes (36) into the r.h.s. and calculates integral with the help of 
(37). Contracting both sides of this relation with the coefficients of the BGW expansion, one finds 

ww = E w = / ^> E 4o» w exp {' s (54) 

R nxn R \ k / 

and the sum in the r.h.s. is the BGW partition function again, only of a different size: 

Z B GWc{n\t) = j y^j-j^ Zbgwc (n + 6 exp ( /? ^ kt k s k j (55) 



Zi°\t) = (52) 
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This relation can be understood as a Fourier transform for the BGW /3-ensemble. It is important 
for two reasons. 

First, it allows one to write t-derivatives of the partition function (correlators) as integrals: say, 
Z BGWc (n\t)\ t=Q = j ^^-Z BGWc (n + 5\s)(^U^)...(^U k -) (56) 



dt kl . . . dtk„ 



whereas without the Fourier transform one could only write 



d" 



dt 



ki 



dt k 



-Z BGWc (n\t) 



d" 



t=o dt kl ...dt km J Vol /3 (n) 



[dU]p a{trU++u us<) 



(57) 



which is largely a symbolical notion: the integrand in the r.h.s. does not depend on the i-variables, 
only the integral does. Thus, the Fourier transform appears to be a convenient tool in this case. 

Second, the Fourier transform allows one to convert the decomposition operator into an integral. 
Indeed, we have 



f 2^A 4k d 2 



ZBGWc(n+\t + )Z B FWc(n-\t ) 



t=0 



(58) 



Substituting each of the BGW partition functions by its Fourier transform, one then finds 

^ - f I KM z — KI*) det i 1 - A ' u+ °T < 59 » 

where m± = n± + ((3 — l)//3 are the shifted sizes. Such a shift looks like a natural property of the 
Fourier transforms for the /3-ensembles. For ,5 = 1, the shift vanishes. 



6 Conclusion 

In this paper, we have obtained two BGW- model-representations of the Selberg models (/3-ensembles), 
elementary constitutients of the conformal blocks, in the pure gauge limit. The spherical Selberg 
model after the PGL turned out to become the character phase of the BGW model, while the elliptic 
Selberg model after the PGL is converted into some other model X. In this paper we succeeded 
in rewriting X as a double BGW model, but we still expect some simpler representations for this 
model. A possible candidate on the role of such a simpler representation is the double-cut BGW 
model, suggested long ago by Dijkgraaf and Vafa [65]. If this is correct, then the double BGW model 
of the present paper is yet another integral representation for this double-cut model. 
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Appendix 1. Properties of Jack polynomials 



The Jack polynomials form an important class of symmetric polynomials, which are often useful in 
calculations with arbitrary /^-ensembles, not only of the Dotsenko-Fateev type. For convenience, in 
this Appendix we list several basic formulas related to the Jack polynomials, which are well-known 
but scattered in the literature. 

Symmetric polynomials 

Symmetric polynomials f(z±, . . . , zn) of given degree degf form a linear space of finite dimension, 

with the basis vectors labeled by the Young diagrams Y = Y\ > Y 2 > Frequently used bases are 

the following: the power sums sy = n^ 5 ^; where Sk{z) = J2i z i are the Newton power sums; the 
elementary symmetric polynomials ey = Y\i e Yi, where e& = the coefficient of x k in Yli(^ + xz i)'i the 
monomial functions my = symmetrization of Ili ^ • The transition matrices between these bases 
have the form 

ei = mi, en = m 2 + 2mn, e 2 = ran, e m = m 3 + 3m 2 i + 6mm, e 2 i = m 2 i + 3mm, e 3 = m m , 
mi = .si, mn = sn/2 - s 2 /2, rn 2 = s 2 , m m = s m /6 - s 2 i/2 + s 3 /3, m 2 i = s 2 i - s 3 , m 3 = s 3 , 
«i = ei, sn = en, s 2 = e n - 2e 2 , s m = e m , s 2i = e m - 2e 2i , s 3 = e m - 3e 2i + 3e 3 , ... 

Often, instead of power sums Sfc(z) = Yli z\ the time-variables tj~ = su/k are used. 
Jack polynomials 

The Jack polynomials Jy are the polynomial eigenfunctions 

WJ Y = J2 Yi ( Yi ~ 2 # " X ) Jy ( 6 °) 

i 

of the VK-like operator 
VF = ^(^(/c-l)-/3/c(/c + l)Js fc — + 2 (fcmg fc+m + / g(fc + m)g fc g 

fc=l Sk fe,m=l ^ Sfc Sm 



<9s 



k+m 



(61) 



normalized with a condition Jy = my + . . . (i.e. the coefficient in front of my is equal to unity). 
Explicitly, several first Jack polynomials with this normalization convention have the form 

Ji(sk) = si 

J 2(Sfc) = S ^J^Y L i Jll(Sk) = \{sl-S 2 ) 

2s 3 + 3f3 Sl s 2 + p 2 s\ (l-P)sis 2 -s 3 + f5s\ 1 3 1 1 

J3(Sfc) " (0 + 1)09 + 2) ' J2l(Sfc) = (/3 + l)(/3 + 2) ' Jlll(sfe) = 6 Sl " 2 Sl * 2 + 3 S3 
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Normalized Jack polynomials: /3-characters 

Jack polynomials are orthogonal with respect to the following scalar product: 



Ja 



Jb)=5ab\\Ja\\ 2 , {sas b ) = \[ 



B j d 



3 3 i 



(62) 



p=0 



which is often called intersection product. The expression for the norm is rather complicated: 



\\Ja\? = ^ 



where Py and Qy are the two "Young normalization factors" 



(63) 



Py= [J (P(Y!-i) + (Yi-j) + p) 



Qy= n (P(Y!-i) + (Yi-j) + l) 



(64) 
(65) 



which have a form of products over the cells (i, j) of the Young diagrams. In applications of the Jack 
polynomials to matrix models, it is often convenient to normalize them with respect to their norm: 



JY 



Jy_ 



(66) 



Such normalized polynomials jy satisfy 



JA 



Jb ) = Sab 



(67) 



and, in general, formulas in terms of jy are simpler than formulas in terms of Jy. However, the 
polynomials jy themselves are more complicated and depend on /3 irrationally. Explicitly, several 
first Jack polynomials with this normalization convention have the form 



jl(Sfc) = VPs!, j 2 ( 



Sk) 



(3{p + 1) s 2 + psj 
2 /3 + 1 



, jll(Sfc) 



/ (3 2 
2j3 + 2 



si - s 2 



>P(P + l){(3 + 2) 2s 3 + 3(3 Sl s 2 + (3 2 s\ 
6 (/3 + l)(/3 + 2) 



j2l(Sfc) 



(2/3 + 1)(3 2 (1 - p) Sl s 2 - s 3 + Psui 



13 + 2 



(0 + 1)03 + 2) 



im(sfe) 



6/? 3 



(2/3 + 1)03 + 



( l 3 1 1 \ 

T)U Sl -2 SlS2 + 3 S3 J 



13 



Particular values of Jack polynomials d R and D R 

Frequently encountered quantities are values of the Jack polynomials at particular points Sk = 
&k,i = (1)0,0, . . .)k and Sk = n (equivalently, tk = 5k,i and tk = n/k). The Jack polynomials have 
the following values at these points: 

J Y {sk = 5 k ,i) = -5— (68) 

Jy(sk = n)= [ ^ (69) 
where [...]y is the so-called " Young- Pochhammer symbol" 

[x]y= II (x-/3(i-l) + (i-l)) (70) 

where Y is any Young diagram and Y 1 stands for its transposed diagram. Note that the Young- 
Pochhammer symbol is a generalization of the classical Pochhammer falling and rising factorials 

(x) n = x(x-l)...(x-n + l) = r( ^ + | ) 1) (71) 

(x)W = x (x + l)...(x + n-l) = + n) (72) 

1 (x) 

which simply correspond to the row and coloumn diagrams Y, respectively. For the normalized Jack 
polynomials jy, the values under consideration are usually denoted through dy and Dy, respectively: 

0\Y\ 

dy = jy{s k = 5 k ,i) = (73) 
Dy{n) = Jy( Sk = n) = ^^= (74) 

This implies, that 

D Y (n) = dy [ ^- (75) 
which is useful to estimate the large-n behaviour of Dy(n). 



14 



Relation to /3-deformed unitary integrals. 

The Jack polynomials are also orthogonal with respect to the integral product 

SASB ) = rTwuT] = 5f if Z (76) 

— 7r — 7r a<b 

which is the /3-deformed unitary integral. The norm is 



Ja 



Jb) = I ^Lj a [u]j b [u+] = s AB 9a V»U (77) 

In J Vol /3 (n) Pa [pn - P + 1]a 



where Vol^rc) = f [dU]p, and Jy[£/] denote the Jack polynomials of power sums Sk = Yla^a °f 
eigenvalues U a = e t( ^ a . In terms of the normalized Jack polynomials, one has 

nxn 

where <5 = ^—^ — is an often encountered shift. 
The /3-deformed Cauchy identity. 

A useful formula is the Cauchy-Stanley identity, sometimes also called the completeness condition: 

exp (X^ SfcSfc ) = 12^ Jr ^ Jr ^ = ^2jd s )jR(s) (79) 
Just as many other formulas, this expansion looks simpler in terms of j-polynomials. 
The /3-deformed IZ integral. 

The character expansion of the /3-deformed IZ integral has a form 

f i du h [ tr (XUYU+) _ST p R Jr[X]J r [Y] _ ^ d R 

J VoW/ e \Q-n [Ma -\D^) 3R[X]3R[Y] m 

nxn n n 

The only origin of n-dependence lies in the explicit factors Dji(n) (equivalently, [j3n] B ). 

The Carlsson-Okounkov shift identity. 

The Carlsson-Okounkov identity has a form 



Ja(sh ~ (~l) fc /3 j - )\j B (s h ~ (~l) k j) 
= II + + + II (-m + p(B' j -i) + (A i -j) + ^ (81) 

(i,j)eA (i,j)eB 

The r.h.s. of this identity has a typical structure of denominators of Nekrasov functions. 
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5d deformations. 



All formulas in this Appendix possess direct generalization to the McDonald polynomials, which 
depend on two parameters. Generalizations of the IZ and BGW integrals can also be defined in this 
way. The McDonald case complements f3 by the "quantum-group" q-deformation and, therefore, is 
relevant for the AGT relation between 5d SYM theories and the conformal blocks of the "quantum" 
Virasoro algebra. For some recent work in this direction, see [16]. In principle, even further defor- 
mations and generalizations may exist, for example, to the Askey- Wilson polynomials [71]. The role 
of Askey- Wilson polynomials in the context of AGT relation remains to be understood. 

Appendix 2. Derivation of Ward identities for the Selberg model 

The Ward identities can be most simply derived by requiring the vanishing of the integral of full 
derivative: 



r N 8 ( N \ 

d Zl ... dz N £ jr-z^ H(z t - Zj f? n ~ L T 'fa - 
1 a=l° Za \i<3 i=l / 



0, s k = J2 z ii k m >0 

i 

(82) 



where L is an auxiliary parameter introduced for the purpose of derivation of the Ward identities 
(we put L = 1 in the end). Whenever possible, we substitute the zero indices of the correlators with 
appropriate powers of N, in order not to work with the to variable. Differentiating the expression 
in the brackets, and using for k m > the relations 

N ft / km 

E z * m+1 £- IE* - = p[-<*m + !)**-. + E s p s ^-p I - ^ 

a=l a i<j \ p=0 J i<j 

N 8 

E ^ dT a Tl Z >= USk - II Z i ( 84 ) 
a=l i i 

N ft N 7 k m + l ( k m \ 

E z » m+1 £- ~ L ) v = v Y. f^i IE* - L ) v = « E Lkm ~ h ^ - Lkm+1 ^ IE* - lt 

a=l a i a=l a i V h=0 / i 

(85) 

N ~ 

Y,z k a ™ +1 — Sl = ls km+l (86) 
„ 1 u a 

a=l 

one finds 

m—l k m 

(u + (k m + 1)(1 - pfjC kl ...k m + ^2 hCk 1 ...k i +k m ...k m - 1 + P^2C kl ...k m _ ukm -p,p + 

i=l p=0 



+ vJ2L km - h C kl ..., km _ lh - L k ™(Ld L )C kl ... km _ x = (87) 

h=0 
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where 

L TV 

C kl ...k m (N) = j d Zl ... dz N H( Zl - z 3 ) 2 ? f[ zf( Zi - L) v s kl ... s km (88) 

are the Selberg integrals with an additional parameter L. Actually, from now on this parameter is 
not needed: it was only useful in derivation of the Ward identities. Using the obvious homogeneity 
property 



(m— 1 \ 
n + j2 k i + uN + vN + n ( n - i )A c kl . 



(89) 



and putting L = 1, one finds that the original integrals Ck 1 ... km (N) = C kl ... km {N) satisfy the 

L = 1 

Ward identities 

m— 1 k m — l 

(u + v + 2/3N+(k m + l)(l- /3)^C kl ... km + J2 hC kl ... ki+km ... km _ 1 +P Yl C ki-k m -i,k m - P , P + 

i=i p=i 



fcm-l / m— 1 \ 

+ v C kl ... km _ uh - [N + Y, h + uN + N(N - l)/3 C k ,... km _, = (90) 

h=l V i=l / 

Note that we explicitly moved the two contributions f3NC kl ,,, km (which arise at particular values 
of p = and p = k m in the third term) into the first term. Similarly, we explicitly moved the 
contributions vC kl _,_ km and vNC kl ,,, km _ 1 into the first and the last terms, respectively. All these 
trivial transformations are necessary to get rid of the zero indices in the correlators and, hence, of 
presence of the to variable in the partition function. 
Because of the obvious formula 



(91) 



the same relations can be rewritten as differential equations known as generalized Virasoro con- 
straints: 



a+b= k h=l 



Zs(h,t 2 , ■■■) = 



m—l 



= f 3 y N+ Yj ki + uN + N(N-l)(3j Z s (h,t 2 ,...), k>0 (92) 

This completes the derivation of the Virasoro constraints for the Selberg model. 

The trick (82) with insertion of a new dimensional parameter L does not work for the elliptic 
Selberg integral, as it does not work in the original Dotsenko-Fateev integrals: dimensionless param- 
eters are present in the both cases. Still analogues of the Virasoro constraints in these both cases 
exist, they will be considered and analyzed elsewhere. 
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